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THE 1973-1984 SOLA! MODULATIOH 07 COSMIC EAT HUCLEI 

* 

M. Garcia-Munox, K.K. Pyle and J.A. Simpson 
Enrico Fermi Institute, University of Chicago 
Chicago, IL, USA 

4222RACT. As a continuation of our program of solar modulation 
studies we have carried out new measurements, with the University of 
Chicago cosmic ray telescope on the Earth satellite IMP-8, of the 
intensity time variations and the energy spectra of galactic cosmic 
ray protons, helium, carbon and oxygen from 1980 through 1984 in- 
cluding the recent solar maximum. In order to test the applicability 
of a steady state model of solar modulation during a period which 
includes times of rapidly changing modulation, we have compared 
these fluxes with the predictions of a conventional model of solar 
modulation which assumes equilibrium between modulation mechanisms • 

It is found that for a reasonable range of variations of the diffu- 
sion coefficient the model predictions can be made to agree with the 
measurements at essentially all times during the studied period. The 
model can account also for the observed hysteresis effects between 
cosmic rays of different rigidities. 

1, Introduction . It has been observed that during the recent solar 
maximum the decreases and increases in cosmic ray intensity due to 
modulation level variations propagate outward from the Sun at ap- 
proximately the solar wind velocity. (McDonald et al., 1981; McKibben, 
Pyle, and Simpson, 1982, 1985; Fillius and Axford, 1985). This fact 
together with the evidence that the radius of the modulation region is 
certainly greater than 30 A.U. and possibly as large as 50 to 150 A.U. 
implies that it takes to the disturbances responsible for these varia- 
tions a time of the order of one year or more for propagation to the 
outer limits of the heliosphere. Therefore, at times of rapidly changing 
modulation, the modulation mechanisms may not be in equilibrium. 

Xhis work is a study of the modulation of cosmic ray nuclei in a 
time interval from 1973 through 1984, including the recent solar maximum. 
Measurements of the intensity time variations and differential energy 
spectra during this period are interpreted in the frame of a conventional 
model of solar modulation which assumes steady state (i.e. equilibrium) 
and absence of drifts. By using empirical diffusion coefficients and 
their time variations it ia found that the model can reproduce at 1 A.U. 
the main features of time intensity variations and differential energy 
spectra of cosmic ray nuclei through the entire solar modulation cycle. 

2. The Modulation Model. The model is described in Evenson et. al.. 
(1983). In it the solar wind velocity, Y, is constant, and equal to 400 
km/sec, and the diffusion coefficient is given by 


k - g P° k(r) ; k(r) - k Q exp [(r-l)/29] (1) 
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where P is par- 
ticle magnet ic 
rigidity, r dis- 
tance from the Sun 
in A ,U. , p par- 
ticle velocity in 
units of velocity 
of light and k 
and o< are func- 
tions of rigidity 
and time. It is 
assumed that the 
local interstellar 
electron differen- 
tial spectrum is 
as deduced from 



synchrotron radiation by Cummings et. aj. . (1973). Comparison with 
measured differential nuclear and electron spectra at 1 A.U. during solar 
minimum leads to the local interstellar nuclear spectra. In this model 
the depth of modulation at a point r is given, for a heliosphere of 
radius R, by the force-field parameter, or modulation parameter 


<£(r) 



( 2 ) 


and thus it is independent of the specific functional form of k(r) 
Therefore, other k(r) functions different from the one selected in (1) 
could lead also to the same modulation depth . 

mn H„i Pe / k ° ( J 983) haVC devel °Ped a time-dependent model of solar 

modulation in which the parameter varying with time is the- frequency of 

outward propagating depressions in k(r) which simulate zones of increased 
scattering associated with solar flare shocks. In our calculations of 

^• ra r- lnt T^ ieS at 1 ^ U » k(r) ®PP ears only in integral form in the 
determination of the modulation depth 
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from 1973 through 1984 
of the 70-95 MeV/n^ 
^uiet time cosmic ray * 
He flux measured at ? 
1A .U . , averaged 
solar rotation 
intervals. The data ^ 
are from the University? 
of Chicago instrument** 
onboard the IMP 8<* 
spacecraft. Note ^hat £ 


the 70-95 MeV/n 4 He £ i 0 4- 


f lux does not contain £ 
anomalous helium, even ^ 
at times of solar 3 
minimum. The anomalous ^ 
helium, when present, 
extends below about 60 
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MeV/n. We have continued measuring the 
differential energy spectra of 10-95 MeV/n 
protons and helium from 1980 through 1984, 
adding to our previous measurements which 
extended up to 1979 (Evenson et. al.. , 1983).^ 
Figure 2 shows the 1977 differential energy • 
spectra at solar minimum, as reported by<^ 
Evenson et aj., (1983), and the 1981 spectra J 
at solar maximum, where we have added the o 
University of New Hampshire data (Webber et -s 
al . . 1983). The curves in the figure are £ 
the fluxes predicted by the modulation 5 
model . Spectral measurements and theoreti- ^ 
cal fits have been obtained in this way * 
from 1977 through 1984. The corresponding 
diffusion coefficients from solar minimum 
1977 through solar maximum 1981, as the 
adjustable parameters, are shown in Figure 


Rigidity Dependence at 1 AU of the Diffusion 
Coefficient from 1977 Through 1981 



P(MV) 


Figure 3 


3. Each of these diffusion coefficients corresponds to a modulation 
depth <f) and gives a modulated helium spectrum with a value for the 70- 
95 MeV/n helium flux. Therefore, there is a one-to-one relationship 
between <f> and the 70-95 MeV/n helium flux. 

The next step is to 
use this relationship and 
use as input the 70-95 
MeV/n He flux, from 1973 
through 1984, in order to 
deduce the model predic- 
tions of the fluxes of 
protons, carbon, and 
oxygen cosmic rays in 
other energy (rigidity) 
intervals through the same 
period and compare these 
predictions with the IMP 8 ^ 

measurements. Figure 4 
shows, for the 1973-1984 «j 
time interval, the com- £ 
parison of model « 
predictions (histograms ) £ 

and measurements (data 2 
points) for 63-95 MeV/n an 
protons, 45-178 MeV/n 
carbon and 53-211 MeV/n 
oxygen. It can be seen 
that in general there is 
good agreement between 
model and measurements. 

The larger disagreements 
can be immediately traced 
to the very well known 

phenomenon of hysteresis Year moiooba 



(Cooper and Simpson, 1979) 
by which higher rigidities 
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recover faster that , 
lower rigidities as can? 
by seen clearly in > 

Figure 4 during the* 
recoveries of the years “ i 
1982 and 1984. However, to 
hys teresis can be c!> 
easily incorporated in w 
our simplified model by'Sj 10 
raising the high^ 
rigidity part of the 3 
diffusion coef f ic ient fi 
relative to the low< 
rigidities during the ^ 
recovery phase. In 3 
this way a good fit to ^ 10 
the 1983 and 1984 
spectra is obtained, as 
is shown, for the particular case of protons and helium, in Figure 5. 

4. Conclusions . Empirical diffusion coefficients can be obtained by 
assuming at all times steady state in the heliosphere and fitting simul- 
taneously the differential energy spectra of protons and helium at 1 AU 
in successive years from solar minimum 1977 through solar maximum 1981. 
These diffusion coefficients define a simple relationship between 70-95 
MeV/n helium flux and depth of modulation. This relationship can be used 
to predict the main features at 1 AU of the intensity time dependence of 
other cosmic ray nucleons in different rigidity intervals from 1973 
through 1984. The hysteresis effects can be accounted for by modifica- 
tion of the diffusion coefficients in selected rigidity intervals. 
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